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On May 25, 2008, the Mars Phoenix Lander successfully landed on the northern arctic 
plains of Mars. An overview of a preliminary reconstruction analysis performed on each en- 
try, descent, and landing phase to assess the performance of Phoenix as it descended is pre- 
sented and a comparison to pre-entry predictions is provided. The landing occurred 21 km 
further downrange than the predicted landing location. Analysis of the flight data revealed 
that the primary cause of Phoenix’s downrange landing was a higher trim total angle of at- 
tack during the hypersonic phase of the entry, which resulted in Phoenix flying a slightly lift- 
ing trajectory. The cause of this higher trim attitude is not known at this time. Parachute 
deployment was 6.4 s later than prediction. This later deployment time was within the varia- 
tions expected and is consistent with a lifting trajectory. The parachute deployment and in- 
flation process occurred as expected with no anomalies identified. The subsequent parachute 
descent and powered terminal landing also behaved as expected. A preliminary reconstruc- 
tion of the landing day atmospheric density profile was found to be lower than the best apri- 
ori prediction, ranging from a few percent less to a maximum of 8%. A comparison of the 
flight reconstructed trajectory parameters shows that the actual Phoenix entry, descent, and 
landing was close to pre-entry predictions. This reconstruction investigation is currently on- 
going and the results to date are in the process of being refined. 


Nomenclature 


a 

= acceleration 

MRO 

= Mars Reconnaissance Orbiter 

C D 

= drag coefficient 

q 

= dynamic pressure 

E 

= atmospheric entry time 

r 

= radius from center of Mars 

EDL 

= entry, descent, and landing 

RSS 

= root sum square 

F 

= force 

S 

= reference area 

IMU 

= Inertial Measurement Unit 

V 

= velocity 

L 

= landing time 

V v 

= vertical velocity 

m 

= mass 

V v 

= vertical velocity 

M 

= Mach number 

Y 

= inertial flight-path angle 

MCS 

MOLA 

= Mars Climate Sounder 
= Mars Orbiter Laser Altimeter 

P 

= atmospheric density 
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I. Introduction 


T HE Mars Phoenix Lander successfully landed on Mars on May 25, 2008 in the northern arctic plains in Vastitas 
Borealis. Its science objective is to search for habitable regions on Mars, continuing the search for the evidence 
of liquid water. Reference 1 provides an overview of the mission. Upon Mars arrival, the capsule containing the 
lander was separated from the cruise stage 7 min prior to entry interface, initiating the entry, descent and landing 
(EDL) phase of the mission. The lander decelerated with the aid of an aeroshell, a supersonic parachute, and 
retrorockets prior to landing safely on the surface. The Phoenix EDL phases and the nominal sequence of events for 
traversing through the atmosphere for the soft landing are illustrated in Figure 1. The Phoenix entry was divided into 
four phases: 1) Entry Preparation - from cruise-stage separation to entry interface, 2) Hypersonic - from entry 
interface to parachute deployment, 3) Parachute - from parachute deployment to lander separation, and 4) Terminal 
Descent - from lander separation to touchdown. Reference 2 provides an overview of the evolution of this Phoenix 
EDL architecture. 


• Final EDL Parameter Update: E-3 hr; Entry State Initialization: E-10 min 

• Cruise Stage Separation: E-7 min 

• Entry Turn Starts: E-6.5 min, Turn completed by E-5 min 

• Entry: E-0 s, L-440 s, 125 km*, r = 3522.2 km, 5.6 km/s, y = -13.0 deg 

• Peak Heating: 46 W/cm 2 Peak Deceleration: 9.2 Earth G 



Entry Prep Phase 


Hypersonic Phase 


Parachute Deployment: E+221 s, L-219 s,12.9 km Mach 1.65 
Heatshield Jettison: E+236 s, L-203 s,11.3 km, 115 m/s 
Leg Deployments: E+246 s, L-191 s 
Radar Activated: E+297 s, L- 143 s 


Parachute Phase 


Landing at 
-4.1 km elevation 
(MOLA relative) 


Lander Separation: E+398 s, L-42 s, 0.94 km, 56 m/s 

• Gravity Turn Start: E+401 s, L-39 s, 0.76 km Terminal Descent 


Constant Velocity Start: E+422 s, L-18 s, 0.051 km, 8 m/s, Phase 
Touchdown: E+440 s, L-0 s, 0 km, Vv=2.0 m/s, Vh=0.7 m/s 
Dust Settling: L+0 to L+15 min 



1 


* Entry altitude referenced to equatorial radius. 
All other altitudes referenced to ground level 



Begin Gyro-Compassing: L+5 min 
• Solar Array Deploy: L+16min 


Figure 1. Phoenix EDL sequence of events. 


The Mars approach maneuvers and targeting procedures were nominal and placed the capsule on the desired 
flight path required for a successful entry. 3,4 The nominal design values at entry interface (radius of 3522.2 km cor- 
responding to a 125 km radial altitude) for the inertial entry velocity and inertial flight-path angle for Phoenix were 
5.6 km/s and -13.0°, respectively. At entry (E) minus 6.5 min, the 3-axis controlled capsule performed a turn-to- 
entry maneuver to attain the desired nominal zero degree angle of attack attitude at entry interface (E-0 s). Early in 
the entry after sensing 0.03 Earth G, the control system deadbands were widened sufficiently to essentially eliminate 
thruster firings due to concerns of a control system instability (i.e., potential reversal) arising from an interaction be- 
tween the aerodynamic flowfield around the capsule and the reaction control system. Reference 5 provides a detailed 
description of this issue. The rationale for the control system deadband settings during the hypersonic phase is pro- 
vided is Ref. 2. As such, during the hypersonic phase, the capsule relied on its inherent aerodynamic stability to 
traverse all the flight regimes (rarefied, transitional, hypersonic, supersonic) until parachute deployment. 6 Although 
Phoenix was designed for a 3 -axis controlled entry, it was actually configured to enter uncontrolled. 2 Consequently, 
Phoenix became the first ballistic non-spinning spacecraft without hypersonic thruster control to land on Mars. 

During the hypersonic phase, the predicted nominal peak deceleration was estimated to be 9.3 Earth G. Para- 
chute deployment was triggered on accelerometer measurements from an Inertial Measurement Unit (IMU) and was 
nominally targeted for an acceleration value of 7.42 m/s 2 (0.756 Earth G), which corresponded to a dynamic pres- 
sure of 490 N/m 2 and a Mach number of 1.65. Heatshield separation was timed to occur 15 s after parachute de- 

2 

American Institute of Aeronautics and Astronautics 

092407 


ployment. Lander leg deployment was timed to occur 10 s after heatshield separation. The radar was timed to be ac- 
tivated 15s after lander leg deployment. Lander separation from the backshell/parachute system was based on radar 
measurements and was nominally to occur at an altitude of 940 m above ground level (AGL). Throttle up and start 
of the gravity turn was timed to initiate 3 s after lander separation, which transitions into a constant velocity phase at 
a vertical velocity of 8 m/s leading to touchdown. 

This paper provides an overview of the preliminary reconstruction analysis performed to date for the Phoenix en- 
try to ascertain how well the actual EDL events compared to the nominal pre-entry sequence of events illustrated in 
Fig. 1 . Accelerometer and gyro data from the IMU and radar data were analyzed to determine the actual event times 
and the as flown trajectory. Each entry phase was examined to assess the performance of Phoenix as it descended to 
quantify how well the pre- entry predictions compared to the actual flight. Lastly, a preliminary comparison of the 
atmospheric density encountered during the entry to the pre-entry profile is presented. 


II. Landing Location 

The navigation accuracy obtained for Phoenix yielded extremely small state errors upon Mars arrival. 4 The final 
orbit determination solution produced the desired nominal inertial entry flight-path angle of -13.0° with a ±3 -g error 
of ±0.007°, which was well within the ±0.15° requirement. The corresponding inertial velocity at entry interface was 
5.6 km/s having a ±3 -a error of ±0.44 m/s. 

Using this final orbit determination solution, a final 
pre-entry predicted landing location (blue circle) was 
calculated that was approximately 17 km north- 
northwest from the desired target of 68.151° North 
Latitude and 233.975° East Longitude (black x) as 
shown in Fig. 2. Note, the Phoenix flight- track went 
from northwest to southeast. The final pre-entry pre- 
dicted landing location was uptrack of the target, be- 
cause TCM-6 was not performed to move it closer. 

The risks associated with performing a TCM so close 
to entry more than outweighed the potential marginal 
gain from a landing closer to the target, since all EDL 
and landing site requirements were already satisfied. 7,8 
The corresponding final pre-entry predicted 99% foot- 
print ellipse was estimated to be 55.1 km by 19.2 km 
having an approach azimuth of 108° (positive from 
North), which was well within the design requirement 
of 1 10 km by 20 km. The actual Phoenix landing loca- 
tion (red diamond) of 68.219° North Latitude and 
234.248° East Longitude was 4.1 km downtrack and 
5.9 km crosstrack from the target, and 21 km down- 
track and 5.3 km crosstrack from the final pre-entry 
prediction as shown in Fig 2. The actual landing loca- 
tion was just inside the final pre-entry predicted 99% 

footprint ellipse, but again was well inside the design requirement of 1 10 km by 20 km. Also depicted in Fig. 2 is the 
landing site elevation relative to the Mars Orbiter Laser Altimeter (MOLA) reference. The Phoenix landing region 
was approximately -4100 m below the MOLA reference. The reason for Phoenix traveling further downtrack and 
crosstrack (towards the northeast) from the predicted landing location is discussed in the next section. However, the 
arrival navigation state can be ruled out as a possible cause, since the desired arrival state was achieved with very 
high accuracy as previously cited. 



East Longitude, deg 

Figure 2. Phoenix landing location (MOLA elevation in 
meters). 


III. Hypersonic Phase 

For the Phoenix entry, nearly the entire deceleration occurs in the hypersonic phase. As such, its performance 
during this phase will dictate the landing location to a very large extent. Analysis of the accelerometer data revealed 
that Phoenix experienced a peak deceleration of 8.5 Earth G occurring at 123 s after entry interface. In comparison, 
the pre-entry nominal peak deceleration predicted was 9.3 Earth G occurring at 121 s after entry interface. Since the 
actual peak deceleration was lower and later than the pre- entry nominal prediction, this situation was an early indi- 
cation that Phoenix’s actual flight was shallower than the prediction, and thus, would lead to a landing further down- 
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track. This shallower flight is observed in Fig. 3, 
which shows Phoenix’s reconstructed relative flight- 
path angle during the hypersonic phase. The flight- 
path angle starts to deviate from the pre-entry predic- 
tion as peak deceleration approaches and remains 
shallower for the duration of the hypersonic phase. 

The reason for this shallower flight was revealed 
by Phoenix’s reconstructed hypersonic attitude his- 
tory. Figure 4 shows the reconstructed total angle of 
attack of Phoenix from entry interface to just prior to 
parachute deployment. As seen, Phoenix trimmed at 
a larger total angle of attack than the pre-entry nomi- 
nal prediction. Early in the entry (from 0 to 75 s), the 
attitude was very similar to the pre-entry nominal 
prediction. However, as the descent continued, 
Phoenix trimmed at a larger total angle of attack. 
From 75 s to 125 s, the trim total angle of attack av- 
eraged approximately 3°, while the estimate for the 
pre-entry nominal prediction was roughly 0.5°. Note, 
although, Phoenix was a ballistic entry, any slight 
radial center of mass offset will result in a small trim 
angle of attack prediction as seen in Fig. 4. From 125 
s to 190 s, the trim total angle of attack was ap- 
proximately 0.9° compared to the pre-entry nominal 
prediction of 0.25°. The exact reason for Phoenix’s 
larger trim total angle of attack is not know at this 
time, and is currently under investigation. However, 
a number of sources that could produce this higher 
trim angle are: 1) a larger radial center of gravity off- 
set than estimated, 2) a mis-prediction of the aerody- 
namics of the capsule, 3) a deformation of the cap- 
sule forebody, or 4) asymmetries in the aeroshell 
forebody. 

The two attitude peaks occurring at approxi- 
mately 77 s and 123 s that increase the total angle of 
attack to approximately 3° are known aerodynamic 
static instabilities which were predicted and expected 
for a blunt body capsule such as Phoenix. Reference 
6 provides an in-depth description of the Phoenix 
capsule aerodynamics and an explanation of the 
cause of these two static instabilities. As the descent 
continue into the supersonic regime (Mach 3 and be- 
low) approaching parachute deployment, the growth 
in the attitude observed from 190 s to 227 s was also 
predicted and expected, and arises from a supersonic 
dynamic instability phenomena that is present in 
blunt bodies. 4 * 6 The pre-entry nominal prediction 
shows this attitude growth increasing from 0.25° to 
2.3°. The actual Phoenix reconstructed attitude is ob- 
served to grow much higher from 0.9° to 4.75° just 
prior to parachute deployment. Note, however, 
Phoenix’s higher trim total angle of attack and larger 
attitude growth approaching parachute deployment 
were within the dispersions characterized through a 
Monte Carlo analysis as well as the design require- 
ments, and are presented later in Section VII. 



Figure 3. Reconstructed relative flight-path angle. 



Figure 4. Reconstructed total angle of attack. 
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Figure 5. Reconstructed bank angle. 


4 

American Institute of Aeronautics and Astronautics 

092407 


Prior to entry, Phoenix was set to a lift-up 
orientation (Bank angle of 180°). The recon- 
structed Bank angle for Phoenix during hyper- 
sonic flight is shown in Fig. 5. As seen, early in 
the entry, the Bank angle remains at the pre- entry 
setting of 180° (pure lift-up). Approaching peak 
deceleration (time of -115 s), however, the Bank 
angle begins to decrease at a constant Bank rate 
of approximately -0.7 deg/s to 107° just prior to 
parachute deployment. During this time, the cap- 
sule rotated an increasing portion of its lift vector 
out of the pure lift-up plane (Bank angle of 180°) 
and into a side- lift orientation (Bank angle of 
107°), thereby exchanging some downtrack 
movement for some crosstrack movement. Con- 
sequently, Phoenix’s larger trim total angle of at- 
tack produced a lifting trajectory that resulted in 
the landing location being 21 km further down- 
track from the pre-entry landing prediction. This 
lifting trajectory is also consistent with producing 
the reconstructed shallower flight-path angle 
(shown in Fig. 3), as well as the later time and 
lower peak deceleration experienced during the 
entry. The reconstructed Bank angle history ex- 
plains the crosstack landing location of 5.3 km 
towards the northeast. The decreasing Bank angle 
during the entry from 180° to 107° indicated a 
slow rotation of the lift vector from pure lift-up 
(i.e., downtrack direction) to a portion into the 
crosstrack direction that is consistent with a 
movement towards the northeast. 

The measured capsule body rates about the 
roll, pitch, and yaw axes during hypersonic flight 
are shown in Figs. 6-8. The roll rate remains near 
zero until peak deceleration when a slow increase 
is observed before settling to a constant value of 
approximately -0.7 deg/s after 125 s. The pitch 
and yaw rates are comparable to the pre-entry 
nominal prediction through peak deceleration. 
The comparison is less favorable during the su- 
personic dynamic instability region; although, the 
larger difference observed is consistent with the 
larger reconstructed attitude growth seen in Fig. 
4. The frequency of the body rates can be used to 
assess the aerodynamics of the capsule, specifi- 
cally its static stability (i.e., the slopes of the 
pitching moment verses angle of attack and yaw- 
ing moment verses side-slip curves). Due to the 
fairly good comparison of the frequency and 
phase of the body pitch and yaw rates to the pre- 
entry nominal prediction is a preliminary indica- 
tion that the aerodynamics generated for Phoe- 
nix 5 6 (static and dynamic coefficients) captured its 
flight behavior reasonably well. A detailed aero- 
dynamic reconstruction is currently underway. 
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Figure 7. Reconstructed body pitch rate. 
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As mentioned previously, Phoenix was designed for a 3 -axis attitude and rate damped controlled entry, but it was 
actually configured to enter uncontrolled due to the widening of the entry control system deadband settings, which es- 
sentially eliminated any thruster firings. 2 As previously stated, this configuration change was due to concerns of a 
control system instability (i.e., potential reversal) arising from an interaction between the aerodynamic flowfield 
around the capsule and the plumes from thruster firings of the reaction control system. 5 In the supersonic regime, 
where the thrusters were most likely to be utilized due to the aerodynamic dynamic instability that would produce 
large attitude and rate growth, the control system deadbands were widened as follows: pitch and yaw attitude dead- 
bands were widened from 5° to 15°; roll attitude deadband was widened from 1.7° to infinity (i.e., no limit); pitch rate 
deadbands widened from 0.5 deg/s to 20 deg/s; yaw rate deadbands widened from 3 deg/s to 20 deg/s; and roll rate 
deadband widened from 0.03 deg/s to 20 deg/s;. As inidicated in Figs. 6-8, the rates that Phoenix experienced during 
its entry would have triggered a large number of thruster firings in all three axes, if the deadbands were not widened, 
increasing the possibility of a control system instability. However, entering with the widened control system dead- 
band configuration, telemetry data indicated that there were no thruster firings during the entry. As a result, Phoenix’s 
predicted inherent aerodynamic stability was sufficient to traverse the entire hypersonic phase of the entry. Overall, 
the hypersonic phase went as expected with no real anomalies identified. Although, Phoenix trimmed at a higher total 
angle of attack than predicted, it was still within the variations expected. However, an investigation is underway to try 
to determine the cause. 


IV. Parachute Phase 

A. Parachute Deployment 

A preliminary reconstruction of the Phoenix para- 
chute deployment, inflation and opening loads was per- 
formed using the accelerometer measurements from the 
IMU. Analysis of the flight data indicated that mortar 
fire occurred 227.84 s after entry interface triggering 
parachute deployment. This mortar fire event is clearly 
seen in Fig. 9, which shows the measured axial accel- 
eration verses time from entry interface. The parachute 
deployment time was 6.4 s later than the pre-entry pre- 
dicted nominal deployment time of 22 1 .44 s after entry 
interface. This later deployment time is consistent with 
the lifting shallower flight trajectory flown by Phoenix. 

The acceleration value just prior to mortar firing was 
7.39 m/s 2 (0.753 Earth G), which was very close to the 
deployment trigger setting of 7.42 m/s 2 (0.756 Earth G). 
The corresponding reconstructed dynamic pressure and 
Mach number at parachute deployment were estimated 
to be 488.9 N/m 2 and 1.7, respectively, which were very 
close to the target values of 490 N/m 2 and Mach 1.65. 
The instantaneous total angle of attack at the parachute 
deployment time of 227.84 s was reconstructed to be 
4.73°, which was well with the design reqirement of 
10°. It appears that parachute deployment was actually 
triggered at the peak of the attitude oscillation observed 
in Fig. 4, as opposed to a possible lower value of the 
cycle, which ranged between 0.25° and 4.75°. The re- 
constructed height at parachute deployment was calcu- 
lated to be 13.3 km (AGL), which was slightly higher 
compared to 12.9 km (AGL) predicted from the pre- 
entry nominal. This higher parachute deployment height 
is again consistent with the lifting trajectory that Phoe- 
nix experienced on landing day. 

The times for line stretch and first peak load from 
mortar fire were 0.75 s and 1.115 s, respectively. These 
times were slightly shorter than the pre- entry estimates 



Figure 9. Measured axial acceleration during para- 
chute deployment. 



Figure 10. Lander RSS body pitch and yaw rates 
during parachute deployment. 

6 


American Institute of Aeronautics and Astronautics 

092407 


of 0.803 s and 1.234 s, respectively. The peak load estimated from the flight data was 44,800 N, which was well be- 
low the flight load requirement of 59,600 N. After peak load, the flight data indicates possible areal oscillations of 
the parachute (i.e., parachute could be breathing - opening and closing, a phenomena associated with disk-gap-band 
parachutes operating at supersonic conditions) before remaining fully open after 230 s. Below Mach 1.4, areal oscil- 
lations are unlikely. 

Figure 10 shows a comparison of the root sum square (RSS) of the capsule body pitch and yaw rates hanging 
under the parachute between the pre-entry nominal prediction and the measured data. Initially, high rates are present 
as expected due to the parachute deployment process, which then damp down very quickly thereafter. As seen, the 
peak measured RSS rates (approx. 90 deg/s) are approximately twice as large as the pre-entry predicted nominal 
(approx. 45 deg/s). This result is a consequence of the larger total angle of attack of 4.75° at parachute deployment 
compared to the much lower value for the pre-entry prediction of 1.8°. However, these high reconstructed rates in 
Fig. 10 are well within the dispersions characterized through a Monte Carlo analysis, and will be presented later in 
Section VII. Overall, the parachute deployment and inflation process went as expected with no anomalies identified. 


B. Parachute Descent 

A reconstruction of the lander descent history dur- 
ing the parachute phase (i.e., from parachute deploy- 
ment to lander separation) is shown in Figs. 11-14. 
The reconstructed lander total angle of attack, pitch 
angle, RSS body pitch and yaw rates, and roll rate 
while hanging under the parachute are compared to 
the pre-entry nominal predication, which was calcu- 
lated using a multi-body parachute/lander simulation. 
Reference 9 provides a detailed description of the de- 
velopment of this multi-body simulation that was util- 
ized during the design phase in an effort to obtain re- 
alistic two-body dynamics of the Phoenix para- 
chute/lander system. The parachute deployment and 
lander separation events from the reconstruction and 
the pre-entry nominal prediction are identified. 

As seen in Figs. 11-14 and discussed in the previ- 
ous subsection, the reconstructed parachute deploy- 
ment time was 6.4 s later than the pre-entry nominal 
prediction. Similarly, the reconstructed time for lander 
separation was also later by 7.8 s (404.84 s verses a 
prediction of 397.04 s). The multi-body simulation 
does reasonably predict the lander total angle of at- 
tack, pitch angle, and RSS rates in comparison to the 
reconstruction. The amplitude and frequency are very 
similar, although they are slightly out of phase (not 
unexpected do to the different parachute deployment 
times). The predicted and reconstructed lander total 
angle of attack (Fig. 11) varied roughly between 5° 
and 25°. The mean total angle of attack in both the 
prediction and reconstruction is seen to increase slow- 
ing during the parachute descent starting from ap- 
proximately 5° after parachute deployment rising to 
approximately 12° prior to lander separation. Figure 
12 compares the lander pitch angle, where -90° points 
vertically down. After parachute deployment, ap- 
proximately 60 s were required for the trajectory to 
turn near vertical. Thereafter, the lander oscillated 
about a pitch angle of -80° (10° from vertical) under 
the parachute. The reason for this 10° offset is be- 
cause the parachute trims at approx. 9° as discussed in 
parachute/lander system, the parachute trim attitude die 



Time from Entry Interface, s 

Figure 11. Lander total angle of attack during the 
parachute descent phase. 



Time from Entry Interface, s 

Figure 12. Lander pitch angle during the parachute 
descent phase. 


Ref. 9. Since the aerodynamics of the parachute dominate the 
dates the trim attitude of the lander as well. 
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The lander RSS body pitch/yaw rates and roll rates 
are shown in Figs. 13 and 14, respectively, during 
parachute descent. After peaking initially during the 
parachute deployment process (as discussed in the 
previous subsection), the RSS rates damp out very 
quickly as observed in Fig. 13. By the times of heat- 
shield separation (15 s after parachute deployment) 
and lander leg deployment (10 s after heatshield sepa- 
ration), the RSS rates have decreased below 30 deg/s 
and remain so for the duration of the parachute de- 
scent phase. Figure 14 shows the roll rate history dur- 
ing parachute descent. The comparison between the 
prediction and reconstructed roll rate is not as good. A 
roll rate peaking at approximately 8 deg/s is observed 
in the reconstruction while descending on the para- 
chute (resulting in the approx. 3.5 revolutions of the 
lander), while there is only a very small roll rate in the 
prediction. However, this reconstructed higher roll 
rate is within the dispersions characterized through a 
Monte Carlo analysis, and is presented later in Section 
VII. 

In future analyses, an attempt will be undertaken 
to determine if the initial high rates experienced by 
Phoenix during parachute deployment can be repro- 
duced by the multi-body simulation, if it was initial- 
ized with the same reconstructed deployment condi- 
tions (e.g., namely the higher parachute deployment 
total angle of attack). In addition, the larger roll rate 
observed in the reconstruction will be further studied 
to determine its cause. Overall, though, the parachute 
phase went as expected with no anomalies identified. 

V. Terminal Descent Phase 

A reconstruction of the lander descent during the 
terminal descent phase (i.e., from lander separation to 
touchdown) is presented in Figs. 15-18. The recon- 
structed and pre-entry nominal prediction lander pitch 
angle and trajectory profiles are shown in Figs. 15 and 
16. Overall, the comparison is observed to be quite 
good. Shortly after separation (927.4 m verses a pre- 
diction of 982.9 m), the lander performed a tip-up 
maneuver (at approx. 800 m) increasing its pitch angle 
to initiate the gravity turn descent to the surface. The 
reconstructed pitch angle (-75°) was slightly lower 
than the prediction (-68°). Also identified in the fig- 
ure is the constant velocity phase just prior to touch- 
down, which began at 52.1 m in the reconstruction 
compared to 56.8 m in the prediction. The telemetry 
data indicated that the backshell avoidance maneu- 
ver 2 , which would divert the lander further upwind to 
provide sufficient separation between it and the para- 
chute in a low or no wind situation, was not em- 
ployed. 

The reconstructed and the predicted lander RSS 
body pitch and yaw rates are shown in Fig. 17. Again, 



Time from Entry Interface, s 


Figure 13. Lander RSS body pitch and yaw rates 
during the parachute descent phase. 



Time from Entry Interface, s 


Figure 14. Lander roll rate during the parachute 
descent phase. 



Figure 15. Lander pitch angle during terminal descent. 
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the comparison between the two is quite good, where 
the magnitude, frequency, and phase of the rate oscil- 
lations are observed to be very similar. Although, the 
reconstruction potentially suggests that the environ- 
ment experienced by Phoenix during its landing was a 
little more benign than in the prediction. The recon- 
structed and predicted lander roll rates are compared 
in Fig. 18. During this phase of descent, the lander is 
designed to rotate to align its orientation to the proper 
landing azimuth for solar array pointing. As seen, the 
magnitude of the roll rate is very similar, although it 
was in the opposite direction compared to the predic- 
tion. This opposite rotation direction is consistent 
with selecting the shorter rotation arc for the 80° roll 
angle change required. Overall, the terminal descent 
phase went as expected with no anomalies identified. 

VI. Preliminary Atmospheric Density 
Reconstruction 

Three weeks prior to entry, daily weather observa- 
tions began of the landing site using the Mars Recon- 
naissance Orbiter (MRO). These observations were 
made to track movements of potential weather sys- 
tems over the landing region that could affect the 
landing. Temperature measurements from the Mars 
Climate Sounder (MCS) instrument 10 on MRO were 
used to make updates to the atmosphere profile used 
in simulations for calculating the landing location. 7 
Over the course of the three weeks, 14 atmosphere 
updates were generated utilizing the MCS tempera- 
ture data. The recommendation of the Atmosphere 
Team was that the ATM012 profile (using MCS data 
from May 18) would best represent the landing day 
atmosphere, since the incremental differences ob- 
served in subsequent profiles where within the pre- 
diction uncertainty and/or were influenced by tran- 
sient dust activity. As are result, ATM012 was used 
to make the final entry parameter update to Phoenix 
and all landing predictions. 7 

Using the acceleration measurement data from the 
IMU during entry, a preliminary atmospheric density 
reconstruction was performed. The atmospheric den- 
sity was calculated using Eq. (1) 

pV 2 

F = -^-C D S = ma (1) 

where F is force, p is atmospheric density, V is the 
reconstructed velocity, C D is the capsule drag coeffi- 
cient, S is the capsule reference area of 5.5155 m 2 , m 
is entry mass of 572.743 kg, and a is the measured 
acceleration from the IMU. All quantities are known 
or measured except for the capsule C D , which can be 
obtained from the aerodynamic database 6 generated 
for Phoenix. However, the assumption made is that 



Figure 16. Lander trajectory profile during terminal 
descent. 



Figure 17. Lander RSS body pitch and yaw rates during 
terminal descent. 



Figure 18. Lander roll rate during terminal descent. 
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there is no or only a small error in the C D prediction. In 
the absence of an onboard air data system to provide 
direct measurement of dynamic pressure (q = pV 2 /2) to 
separate the density and C D terms in Eq. (1), any errors 
in the pre-entry predicted value of C D do have a direct 
impact on the atmospheric density reconstruction. 

However, historically, the C D for blunt bodies during 
hypersonic flight has been predicted well (within a 
couple percent), since the aerodynamics are entirely 
forebody dominated. 11 The uncertainty in the accelera- 
tion measurement (a) must also be considered in 
evaluating the accuracy of the reconstruction. How- 
ever, the specific Litton-200 IMU used by Phoenix 
was very accurate. 

Given this assumption, the atmospheric density en- 
countered during the Phoenix entry was reconstructed 
and is shown in Fig. 17, where it is plotted as a fraction 
of the ATM012 profile; that is, if the reconstructed 
density were identical to ATM012, a vertical line at a 
value of 1 .00 would result. As seen, between approxi- 
mately 13 km and 48 km, which is where nearly all of the deceleration occurs, the reconstructed density is less than 
ATM012 ranging from a few percent lower to a maximum reduction of 8% at 29 km. This 8% reduction corre- 
sponds to approximately a 1.5 -g low variation from the atmosphere dispersions generated for Phoenix. 12 Note, the 
atmosphere reconstruction stops at parachute deployment (approximately 13 km), since the confidence in the para- 
chute C D value is low resulting in an unreliable density estimate. 

Also shown in Fig. 17 is ATM021, which was generated after landing using MCS temperature data measured 
approximately two hours prior to entry on landing day May 25. This temperature data set would produce the best es- 
timate of the landing day atmosphere profile; however, the time required to downlink and process this temperature 
data precluded its use prior to entry. Between 13 km to 30 km, ATM021 is similar to ATM012. Above 30 km, 
ATM021 is denser than ATM012. As such, even if ATM021 was available prior to entry, the reconstructed density 
would still be less than the best apriori prediction for the landing day atmosphere profile. Shortly after landing, a 
measurement of the surface pressure and temperature was made and recorded. In future analysis, this preliminary 
density reconstruction will be refined using this additional data by propagating the density up from the surface. 

This lower atmospheric density experienced by Phoenix compared to ATM012 is also consistent with landing 
further downrange. To ascertain the sensitivity of this lower density on the landing location, it was substituted for 
ATM012 in the same simulation utilized for predicting the landing day location. The resulting shift in landing loca- 
tion was only 4.2 km further downtrack. Consequently, this reconstructed lower atmospheric density experienced by 
Phoenix was not the root cause and was only a small contributor to the landing being 21 km further downtrack of the 
prediction (see Section II). As a result, the predominant reason for Phoenix’s further downrange landing was due to 
a larger trim total angle of attack (see Section III) resulting in a lifting trajectory that carried it downtrack. 

VII. Reconstruction Performance Overview 

The robustness of the Phoenix entry design was assessed through a Monte Carlo dispersion analysis to statisti- 
cally determine its performance to off-nominal conditions. References 7, 9, and 13 describe the details of the simula- 
tion developed to model the Phoenix EDF phases and the assumptions of the dispersion analysis. In the Monte Carlo 
dispersion analysis, 2000 random off-nominal cases were simulated varying uncertainties in all the EDF parameters. 
For example, variations in the initial state vector, capsule mass properties (mass, center of gravity, inertia), initial at- 
titude and attitude rates, hypersonic aerodynamic coefficients, atmospheric density and winds, parachute aerody- 
namics, descent engine thrust, etc were considered. Results from the final pre-entry prediction using the multi-body 
simulation 9 are shown in Table 1, where entry conditions are given at EDF events during the descent in terms of 
their statistical mean and 99% high and low variations. These results were the best apriori estimates of the expected 
entry conditions and their variations. 

Also listed in Table 1 for comparison to the pre-entry predictions are trajectory parameters that have been pre- 
liminarily reconstructed thus far from the Phoenix flight data. As seen, majority of the reconstructed parameters are 
close to the mean predictions. For example, the actual parachute deployment dynamic pressure and fuel used during 
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Density/Density ATMQ12 


Figure 17. Preliminary atmospheric density 
reconstruction. 


the entry and descent were at the mean predicted value. There are, however, a few reconstructed parameters that ap- 
proach the 99% high or low variations or even exceed them slightly (e. g., peak deceleration, attitude at peak heat 
rate). Although, all of the reconstructed parameters are well within the listed EDL design requirements. 7 Overall, the 
actual Phoenix entry, descent, and landing was near nominal and close to pre-entry predictions. The reconstruction 
investigation is ongoing and the results to date are in the process of being refined. 


Table la: Comparison of Phoenix Pre-Entry Predication and Preliminary Reconstruction Results I 


Parameter 

Units 

99% 

Low 

Mean 

99% 

High 

Design 

Requirement 

Preliminary 

Reconstructed 

Hypersonic Flight 







Peak Heating Rate 

W/cm 2 

44.7 

46.1 

47.4 

<64.0 


Attitude at Peak Heat Rate 

deg 

0.1 

0.6 

1.5 

< 10.0 

2.6 

Peak Deceleration 

Earth G 

8.7 

9.3 

9.8 

< 13.0 

8.5 

Total Heat Load 

J/cm 2 

2305.9 

2405.8 

2509.4 

<3320 


Parachute Deployment 







Time from Entry 

s 

211.3 

219.9 

229.4 


227.8 

Height (AGL) 

km 

10.0 

12.7 

15.7 


13.3 

Wind-Relative Velocity 

m/s 

330.9 

368.3 

419.1 


387.6 

Mach Number 


1.45 

1.64 

1.89 

1.1 <M<2.13 

1.7 

Dynamic Pressure 

N/m 2 

452.5 

490.5 

532.9 

300 < q < 560 

488.9 

Relative Flight-Path Angle 

deg 

-33.9 

-28.3 

-22.6 


-27.1 

Relative Azimuth 

deg 

106.2 

111.7 

116.8 


109.8 

Total Angle of Attack 

deg 

0.2 

2.3 

7.3 

< 10.0 

4.73 

Sensed Acceleration 

Earth G 

0.744 

0.770 

0.792 


0.753 

RSS Pitch/Yaw Rates 

deg/s 

0.4 

4.5 

16.1 


6.9 

Roll Rate 

deg/s 

-1.8 

0.0 

1.7 


-0.83 

Max RSS Pitch/Yaw Rates 







thru parachute deployment 

deg/s 

25.3 

67.9 

171.1 


93.1 

Heatshield Jettison 







Height (AGL) 

km 

8.3 

11.1 

13.9 


11.6 

Wind-Relative Velocity 

m/s 

97.7 

116.3 

141.5 


123.0 

Relative Flight-Path Angle 

deg 

-52.7 

-44.1 

-32.6 


-40.5 

Dynamic Pressure 

N/m 2 

45.9 

56.4 

68.3 



Mach number 


0.43 

0.51 

0.62 

<0.8 


Total Angle of Attack 

deg 

0.7 

7.9 

20.8 


7.8 

Sensed Acceleration 

Earth G 

0.73 

0.83 

0.95 



RSS Pitch/Yaw Rates 

deg/s 

1.3 

15.5 

48.3 

< 100.0 

16.5 

Roll Rate 

deg/s 

-3.8 

-0.1 

3.4 


5.5 
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Table lb: Comparison of Phoenix Pre-Entry Predication and Preliminary Reconstruction Results II 


Parameter 

Units 

99% 

Low 

Mean 

99% 

High 

Design 

Requirement 

Preliminary 

Reconstructed 

Lander Leg Deployment 

Height (AGL) 

km 

7.5 

10.2 

13.1 


10.9 

Wind-Relative Velocity 

m/s 

71.5 

83.3 

99.0 


91.3 

Dynamic Pressure 

N/m 2 

25.8 

31.1 

36.8 



Mach number 


0.31 

0.34 

0.44 



Total Angle of Attack 

deg 

0.7 

7.9 

20.8 


7.0 

RSS Pitch/Yaw Rates 

deg/s 

1.0 

12.6 

40.1 

< 100.0 

5.8 

Roll Rate 

deg/s 

-3.9 

-0.1 

19.6 


6.4 

Lander Separation 

Time from Entry 

s 

348.1 

392.3 

435.3 


404.9 

Time from Parachute Deploy 

s 

134.4 

172.4 

211.5 


177.1 

Height (AGL) 

m 

922.9 

982.9 

1249.4 


925 

Wind-Relative Velocity 

m/s 

44.3 

54.3 

64.4 


56.7 

Relative Flight-Path Angle 

deg 

-86.6 

-78.6 

-64.8 


-75.7 

Mach number 


0.19 

0.22 

0.26 



Total Angle of Attack 

deg 

1.9 

11.9 

25.0 


19 

Pitch Angle 

deg 

-88.5 

-77.5 

-62.5 


-79.2 

RSS Pitch/Yaw Rates 

deg/s 

1.2 

11.8 

32.9 

< 100.0 

4.8 

Roll Rate 

deg/s 

-5.1 

-0.9 

3.4 


-0.97 

Tip-Up 

Height 

m 

892.4 

951.9 

1213.6 


897 

Gravity Turn 

Height 

m 

746.9 

805.8 

1045.1 


720 

Constant Velocity 

Height 

m 

45.9 

51.9 

57.7 


52.1 

Landing 

Time from Entry 

s 

392.8 

436.2 

479.0 


446.1 

Time from Parachute Deploy 

s 

177.4 

216.3 

254.8 


218.3 

Time from Lander Separation 

s 

40.5 

43.9 

48.7 


41.2 

Relative Vertical Velocity 

m/s 

1.4 

2.1 

2.7 

1 .0 < V v < 3 .4 

2.38 

Relative Horizontal Velocity 

m/s 

0.0 

0.5 

1.2 

< 1.4 

0.06 

Fuel Usage During Entry & Descent 

kg 

35.3 

37.6 

41.5 


37.6 

Landing Footprint size 

km 




110x20 

55.1 x 19.2 


VIII. Conclusion 

The Mars Phoenix Lander successfully landed on Mars on May 25, 2008. A reconstruction of Phoenix’s entry 
was performed. The landing was 21 km further downrange with a northern crosstrack of 5.3 km than the predicted 
landing location. Analysis of the flight data revealed that the cause of Phoenix’s downtrack and crosstrack landing 
was a higher trim total angle of attack during the hypersonic phase of the entry, which resulted in Phoenix flying a 
lifting trajectory. Approximately a 3° trim total angle of attack was observed compared to a prediction of 0.25°. The 
cause of this higher trim attitude is not known at this time. Parachute deployment was 6.4 s later than the pre-entry 
prediction. This later deployment time was within the variations expected and is consistent with a lifting trajectory. 
The parachute deployment and inflation process went as expected with no anomalies identified. The peak parachute 
inflation load was estimated to be 44,800 N, which was well below the 59,600 N flight load requirement. During the 
subsequent descent, the oscillations of the lander under the parachute were close to pre-entry predictions. Similarly, 
the powered terminal descent phase behaved as expect. A preliminary reconstruction of the landing day atmospheric 
density profile was found to be lower than the a prior prediction, ranging from a few percent less to a maximum of 
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8%. Overall, the preliminary reconstructed entry, descent, and landing trajectory parameters were found to be close 
to mean pre-entry predictions or within the expected variation. The reconstruction effort is currently ongoing and the 
results to date are in the process of being refined. 

Acknowledgments 

The authors would like to specifically acknowledge the effort of Erik Bailey at the Jet Propulsion Laboratory for 
processing the Phoenix telemetry data into a form that was usable for this reconstruction investigation. Additionally, 
the authors would like to recognize the contributions of the entire Phoenix EDL team at the Jet Propulsion Labora- 
tory (Doug Adams, Erik Bailey, Gene Bonfiglio, Jim Chase, Curtis Chen, Ben Cichy, Dan Eldred, Richard Korn- 
feld, Paul Laufer, Mike Lisano, Rob Manning, Dana Runge, Dara Sabahi, David Skulsky, Cris Windoffer), at Lock- 
heed Martin Space Systems (Tim Gasparrini, Brad Haack, Mark Johnson, Tim Linn, Tim Priser, Jay St. Pierre), and 
at NASA Langley Research Center (Jody Davis, Artem Dyakonov, Chris Glass, Karl Edquist, Ray Mineck, John 
Van Norman) for the overall success of the Phoenix landing. A portion of this research was carried out at the Jet 
Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space 
Administration. Reference herein to any specific commercial product, process, or service by trade name, trademark, 
manufacturer, or otherwise, does not constitute or imply its endorsement by the United States Government or the Jet 
Propulsion Laboratory, California Institute of Technology. 

References 

^hotwell, R. “Phoenix - The First Mars Scout Mission (A Midterm Report),” Acta Astronautica 57, July-October 2005, pp. 
121-134. 

2 Grover, M. R., Cichy, B. D., and Desai, P. N., “Overview of the Phoenix Entry, Descent, and Landing System Architecture,” 
AIAA Paper 2008-7218, Proceedings of AIAA Guidance, Navigation, and Control Conference, Honolulu, HI, August 2008. 

3 Portock, B. M., Kruizinga, G., Bhat, R., et al., “Navigation Challenges of the Mars Phoenix Lander Mission,” AIAA Paper 
2008-7214, Proceedings of AIAA Guidance, Navigation, and Control Conference, Honolulu, HI, August 2008. 

4 Ryne, M. S., Graat, E. J., Lau, E., et al., “Orbit Determination for the 2007 Mars Phoenix Lander Mission,” AIAA Paper 
2008-7215, Proceedings of AIAA Guidance, Navigation, and Control Conference, Honolulu, HI, August 2008. 

5 Dyakonov, A. A., Glass, C. E., Desai, P. N., and Van Norman, J. W., “Aerodynamic Interference Effects Due to Reaction 
Control System for the Phoenix Entry Capsule,” AIAA Paper 2008-7220, Proceedings of AIAA Guidance, Navigation, and Con- 
trol Conference, Honolulu, HI, August 2008. 

6 Edquist, K. T., Desai, P. N., and Schoenenberger, M., “Aerodynamics Prediction for the Mars Phoenix Entry Capsule,” 
AIAA Paper 2008-7219, Proceedings of AIAA Guidance, Navigation, and Control Conference, Honolulu, HI, August 2008. 

7 Prince, J. L., Desai, P. N., Queen, E. M., Grover, M. R., “Entry, Descent, and Landing Operations Analysis for the Mars 
Phoenix Lander” AIAA Paper 2008-7349, Proceedings of AIAA Guidance, Navigation, and Control Conference, Honolulu, HI, 
August 2008. 

8 Bonfiglio, E. L., Spencer, D. A., Arvidson, Stein, T., and Heet, T., “Landing Site Dispersion Analysis and Statistical As- 
sessment for the Mars Phoenix Lander” AIAA Paper 2008-7348, Proceedings of AIAA Guidance, Navigation, and Control Con- 
ference, Honolulu, HI, August 2008. 

9 Queen, E. M., Desai, P. N., Prince, J. L, and Cruz. J. R., “Multi-Body Modeling and Simulation for Mars Phoenix Lander 
Entry, Descent, and Landing,” AIAA Paper 2008-7347, Proceedings of AIAA Guidance, Navigation, and Control Conference, 
Honolulu, HI, August 2008. 

10 McCleese, D. J., Schofield, J. T„ Taylor, F. W„ Calcutt, S. B„ Foote, M. C., Kass, D. M„ Leovy, C. B„ Paige, D. A., Read, 
P. L., Zurek, R. W., “Mars Climate Sounder: An investigation of thermal and water vapor structure, dust and condensate distribu- 
tions in the atmosphere, and energy balance of the polar regions,” Journal of Geophysical Research , Vol. 112(E5), DOI: 
1 0.1 029/2006JE002790, 2007. 

1 Edquist, K. T., “Computations of Viking Lander Capsule Hypersonic Aerodynamics with Comparisons to Ground and 
Flight Data,” AIAA Paper 2006-6137, AIAA Atmospheric Flight Mechanics Conference and Exhibit, , Keystone, CO, August 
2006. 

12 Tamppari, L., J. Bames, E. Bonfiglio, B. A. Cantor, A. J. Friedson, A. Ghosh, M. R. Grover, D. Kass, T. Z. Martin, M. T. 
Mellon, T. Michaels, J. R. Murphy, S. Rafkin, M. D. Smith, G. Tsuyuki, D. Tyler, and M. J. Wolff (2008), “The Atmospheric 
Environment Expected for the Phoenix Landed Season and Location,” Journal of Geophysical Research, 
doi: 1 0. 1 029/2007JE003034, in press. 

13 Prince, J. L., Desai, P. N., Queen, E. M., Grover, M. R., “Mars Phoenix Entry, Descent, and Landing Simulation Design 
and Modeling Analysis” AIAA Paper 2008-7507, Proceedings of AIAA Guidance, Navigation, and Control Conference, Hono- 
lulu, HI, August 2008. 


13 

American Institute of Aeronautics and Astronautics 

092407 



